DIFFERENTIALLY COHERENT COMBINING FOR ELECTRONIC ARTICLE 
SURVEILLANCE SYSTEMS 


CROSS REFERENCES TO RELATED APPLICATIONS 
5 This application claims the benefit of U.S. Provisional Application No. 60/267,886, 

filed February 8, 2001. 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR 
DEVELOPMENT 

p 10 Not Applicable 

§Jl BACKGROUND OF THE INVENTION 

f I Field of the Invention 

fy This invention relates to electronic article surveillance receivers, and more 

L 1 5 particularly, to signal processing and detection techniques for an electronic article 

nj surveillance receiver. 

o 

f i. Description of the Related Art 

O Electronic article surveillance (EAS) systems, such as disclosed in U.S. Patent No. 

4,5 10,489, transmit an electromagnetic signal into an interrogation zone. 

20 Magnetomechanical EAS tags in the interrogation zone respond to the transmitted signal with 
a response signal that is detected by a corresponding EAS receiver. Pulsed 
magnetomechanical EAS systems have receivers, such as ULTRA*MAX receivers sold by 
Sensormatic Electronics Corporation, Boca Raton, Florida, that utilize noncoherent detection 
and a highly nonlinear post detection combining algorithm in processing the received signals. 

25 To improve processing gain, phase information present in the received signal can be utilized 
in detection. 

BRIEF SUMMARY OF THE INVENTION 
A system and method for differential coherent combining of received signals in an 
30 electronic article surveillance receiver is provided. The systems includes receiving a receive 
signal including a first component of an electronic article surveillance tag response and a 
second component of noise. Next the receive signal is filtered with a plurality of filters each 


C4-988 


1 


having a preselected bandwidth and a preselected center frequency. The output of each of 
said plurality of filters are sampled to form a plurality of filtered samples. Each of the 
plurality of filtered samples are combined by diversity averaging. A quadratic detector 
detects each of the plurality of filtered samples by squaring the diversity combined samples 
and summing to arrive at a differentially coherent combined signal. 

The system may further compare the differentially coherent combined signal to a 
preselected threshold and provide an output signal associated with said comparison. The 
output signal may trigger an alarm or other selected reaction. 

Objectives, advantages, and applications of the present invention will be made 
apparent by the following detailed description of embodiments of the invention. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 

Figure 1 is a block diagram of a conventional EAS transmitter. 

Figure 2 is a plot of a transmit signal and tag response signal. 

Figure 3 is a block diagram of a conventional matched filter detector. 

Figure 4 is a block diagram of a conventional quadrature matched filter detector. 

Figure 5 is a block diagram of an implementation of a bank the quadrature matched 
filters shown in Fig. 4. 

Figure 6 is a block diagram of the bank of filters of Fig. 5 with conventional initial 
hit/validation combining. 

Figure 7 is a plot of receiver operating characteristics of coherent and noncoherent 
detection. 

Figure 8 is a block diagram illustrating the inventive detector using differential 
coherent combining. 

Figure 9 is flow chart of the outlier discrimination algorithm. 

DETAILED DESCRIPTION OF THE INVENTION 
Referring to Fig. 1, a conventional pulsed EAS transmitter is illustrated, such as that 
sold under the name ULTRA* MAX by Sensormatic Electronics Corporation. Sequence 
Controller 2 is typically a state machine that executes in software. It is responsible for 
frequency hopping and phase flipping the transmit signal so that tags of various center 
frequencies and physical orientations are adequately excited by the transmitter. The 
frequency control signal f(t) takes on one of three values. When f(t) = 0, then the nominal 
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center frequency, such as 58,000 Hz, is transmitted. When f(t) = 1, then the high frequency, 
such as 58,200 Hz is transmitted. If f(t) = -1, the low frequency, such as 57,800 Hz is 
transmitted. The phase control signal p(t) takes on one of two values, p(t) = 1 or p(t) = -1 . 
This controls the polarity of the transmit antennas 4, either aiding or opposing. The carrier 
signal is typically a phase locked loop based oscillator that includes a voltage controlled 
oscillator 6 that is modulated by the phase and frequency control inputs 8. The carrier signal 
c(t) can be denoted: 

C(t) = p(t)-sin(2-7t- (f c + f(t))-t + 6), 

where 9 is an arbitrary phase angle that depends on the hardware. The carrier signal is 
combined 10 with a baseband pulse train m(t) before being amplified 12 

The receive signal is processed by an analog front end, sampled by an analog to 
digital converter (ADC), and compared to a threshold. The threshold is set by estimating the 
noise floor of the receiver, then determining some suitable signal to noise ratio to achieve a 
good trade off between detection probability, P de t, and false alarm probability, P fa . The 
sequence controller 2 would typically produce frequency and phase control signals as shown 
in Fig. 1. When a signal is initially detected based on the threshold test (known as an "initial 
hit"), the sequence controller 2 "locks" the transmitter phase and frequency values for a 
"validation sequence". The validation sequence is usually around six transmit bursts long. 
During this validation sequence the system basically verifies that the signal continues to be 
above the threshold. 

There are two modes of operation for a magnetomechanical tag, such as an 
ULTRA*MAX tag as disclosed in the '489 patent, linear and nonlinear. For the linear model, 
the tag behaves as a simple second order resonant filter with impulse response: 

h(t) = Ao-e at -sin(2-7*- f„ -t + 9) 

where A 0 is the amplitude of the tag response, f„ is the natural frequency of the tag, and a is 
the exponential damping coefficient of the tag. Fig. 2, shows a plot of a transmit signal 14 
and the tag response signal 16 when the tag operates linearly. 

The nonlinear model is more closely coupled to the mechanics of the tag itself. The tag 
becomes nonlinear when it is overdriven by the transmitter. In this case, the resonator(s) 
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within the cavity vibrate so hard that they begin to bounce off the interior walls of the cavity. 
In this mode, the behavior is analogous to the ball inside the pinball machine. Very small 
changes in initial conditions of the resonator result in large changes in the phase and 
amplitude of the final tag ring down. This is an example of the nonlinear dynamics known as 
5 chaos. Although this nonlinear response will be mentioned briefly, the present invention is 
primarily concerned with detection of the tag when it is in the region of linear behavior. 
Thus, unless specifically called out, the remainder of this description refers to tag response 
that is linear. 

The signal from the receive antenna when a tag is present is the sum of the tag's 

1 0 natural response to the transmit signal plus the additive noise due to the environment. 

ULTRA*MAX systems operating around 60000 Hz preside in a low frequency atmospheric 
noise environment. The statistical characteristics of atmospheric noise in this region is close 
to Gaussian, but somewhat more impulsive (i.e., a symmetric cc-stable distribution with 
characteristic exponent near, but less than, 2.0). 

1 5 In addition to atmospheric noise, the 60000 hertz spectrum is filled with man-made 

noise sources in a typical office/retail environment. These man-made sources are 
predominantly narrowband, and almost always very non-Gaussian. However, when many of 
these sources are combined with no single dominant source, the sum approaches a normal 
distribution (due to the Central Limit Theorem). 

20 The classical assumption of detection in additive white Gaussian noise is used herein. 

The "white" portion of this assumption is reasonable since the receiver input bandwidth of 
3000 to 5000 hertz is much larger than the signal bandwidth. The Gaussian assumption is 
justified as follows. 

Where atmospheric noise dominates, the distribution is known to be close to 

25 Gaussian. Likewise, where there are a large number of independent interference sources the 
distribution is close to Gaussian due to the Central Limit Theorem. If the impulsiveness of 
the low frequency atmospheric noise were taken into account, then the optimum detector 
could be shown to be a matched filter preceded by a memoryless nonlinearity. The optimum 
nonlinearity can be derived using the concept of influence functions. Although this is 

30 generally very untractable, there are several simple nonlinearities that come close to it in 
performance. To design a robust detector we need to include some form of nonlinearity. 
When there is a small number of dominant noise sources we include other filtering to deal 
with these. For example, narrow band jamming is removed by notch filters or a reference 
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based least means square canceller. After these noise sources have been filtered out, the 
remaining noise is close to Gaussian. Although many real installations may deviate from the 
Gaussian model, it provides a controlled, objective set of conditions with which to compare 
various detection techniques. 

Referring to Fig. 3, when the signal of interest is completely known a matched filter is 
the optimum detector. In our case, say we knew the resonant frequency of the tag and its 
precise phase angle when ringing down. The signal we're trying to detect is 

s(t) = A- e at -sin(2-7t- f n t + 0). 

Then the matched filter is simply the time reversed (and delayed for causality) signal, s(T r - 
t) at 1 8. The matched filter output is sampled 20 at the end of the receive window, T r , and 
compared to the threshold 22. A decision signal can be sent depending on the results of the 
comparison to the threshold. The decision can be a signal to sound an alarm or to take some 
other action. Note that we do not have to know the amplitude, A. This is because the 
matched filter is a "uniformly most powerful test" with regard to this parameter. This 
comment applies to all the variations of matched filters discussed below. 

Referring to Fig. 4, when the signal of interest is completely known except for its 
phase 0, then the optimum detector is the quadrature matched filter (QMF). QMF is also 
known as noncoherent detection, since the receiver is not phase coherent with the received 
signal. On the other hand, the matched filter is a coherent detector, since the phase of the 
receiver is coherent with the received signal. The receive signal r(t) which includes noise and 
the desired signal s(t) is filtered by s(T r - 1) at 24 as in the matched filter, and again slightly 
shifted in phase by n/2 at 25. The outputs of 24 and 25 are sampled at 29, squared at 26 and 
27, respectively, combined at 28, and compared to the threshold 30. 

Referring to Fig. 5, when the signal of interest is completely known except for its 
frequency f n and phase 9, then the optimum detector is a bank of quadrature matched filters 
(QMFB). A quadrature matched filter bank can be implemented as a plurality of quadrature 
matched filters 40, 42, and 44, which correlate to quadrature matched filters with center 
frequencies off, f 2 through f m , respectively. The outputs of the quadrature matched filters 
are summed at 46 and compared to a threshold at 48. 

Referring to Fig. 6, often the signal to noise ratio does not allow for the desired 
performance, i.e., low enough false alarm probability P fa with high enough detection 
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probability P det . In this case one form or another of diversity may be available to improve the 
SNR, thereby reaching performance goals. Systems such as ULTRA* MAX use time 
diversity, averaging over multiple receive windows to reduce the effects of noise. The 
textbook method for doing this with a quadrature matched filter bank is to average the QMFB 
output over many receive windows and perform a threshold test. For white Gaussian noise, 
the noise in different receive windows is uncorrelated and therefore its effects can be reduced 
by averaging. Asymptotically, the noise can be reduced 1.5dB for every doubling of the 
number of receive windows averaged. However, using coherent detection 3.0dB of noise 
reduction can be achieved for every doubling of the number of receive windows averaged. 
This is a significant difference and is an important feature of the present invention. 

Present EAS systems using nonlinear post detection combining is illustrated by the 
initial hit/validation diversity combiner 50. The resulting detection statistic is compared to an 
estimate of the noise floor. If a signal to noise ratio criteria is met the system will go into 
validation. At this point the sequence controller 2, shown in Fig. 1, locks to the transmitter 
configuration which passed the initial hit threshold test. The transmitter does a number of 
additional bursts N, typically about six. If all N of the receive samples pass the threshold test, 
then the system alarms. 

This validation sequence is in effect a form of post detection combining, albeit a very 
nonlinear one. It can be referred to it as a "voting" combiner, where a certain percentage of 
the threshold tests must pass, for example, this may require 100% pass, for a unanimous 
decision. 

To analyze the performance of the conventional detection scheme, specifically the 
noncoherent detection with "initial hit/validation" type post detection combining, we assume 
a Neyman-Pearson type criteria, i.e., we choose an acceptable level for the false alarm rate 
P fa , then determine our probability of detection P det verses SNR. Receiver operating 
characteristics for coherent and noncoherent detection, as well known in the art, is shown in 
Fig. 7. 

First, the probability of passing the threshold test on a single receive test statistic when 
in fact there is no tag signal present is denoted as P&, the probability of false validation. A 
validation sequence would follow in which all N test statistics would have to be above the 
threshold. Using the independence of the receive samples we have 

Pfa=Pfv (N+1) - 
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Likewise, P lh is the probability of passing the threshold test when there is in fact a tag signal 
present. Again using independence, the probability of detection is 


Pdet = P lh (N+l) . 

Now, we choose N = 3 and P fa = 10" 8 . Solving, we get Pfv = 10' 2 . Assume that the threshold 
is set for 12dB. Then using the curves in Fig. 7 for noncoherent detection, P lh = 0.992. Then 
calculating P de t = 0.968. 

Notice that if only one receive sample at P fa = 10" 8 and 12dB SNR, then P det = 0.35. 
To achieve P det = 0.968 we would have needed 14.8 dB SNR. This difference, 14.8 dB - 12 
dB = 2.8 dB, represents the processing gain due to the "unanimous vote" combining scheme 
used in the conventional receiver. 

It is apparent that a great deal of information is being lost by ignoring the signal's 
phase. The data is reduced beyond the point of a sufficient statistic (we no longer satisfy the 
sufficiency requirement fundamental to detection theory). The present invention recovers 
this lost information. The result is based on the linearity of the tag model, and transposing 
the order of linear post detection combining and noncoherent detection. 

Since the tag signal is linear, then given a set of initial conditions and parameters a, 
and f n , its response is determined. For any given tag in a given orientation, its parameters are 
fixed. Therefore, if the transmitter function is the same for every transmit burst, then the 
tag's initial conditions when the transmitter shuts off will be the same, and the tag's natural 
response will be the same. That is, the tag signal's amplitude A and phase 0 will be fixed. 

This turns out to be true over short durations of time even when the tag is in motion. 
In other words, when the tag passes through the interrogation zone at one meter per second in 
a set orientation, its phase changes very little. Its amplitude changes relative to the amount ol 
transmitter field it is excited by. However, given that the transmitter repetition rate is about 
90 hertz (one burst every 1 1 milliseconds) the tag can only move 1 1 millimeters in this time. 
Over short periods of time the tag's amplitude is relatively stable. 

The fact that the tag signal's amplitude and phase are approximately equal from one 
receive window to the next is valuable information. The exact value of the signal's phase is 
not known, but we know that the differential of the phase angle is nearly zero. To take 
advantage of this, diversity combining can be implemented in front of the quadrature 
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detector. This takes advantage of the 3.0dB per doubling processing gain of coherent 
combining without actually knowing the signal's phase. 

Note that to accomplish this processing gain, the system must do away with the 
concepts of initial hit and validation. Instead, the sequence controller portion of the 
transmitter must now send N identical transmit bursts in a row prior to any decision being 
made by the detector. This is analogous to the fixed length dwell concept used in radar 
systems. 

Referring to Fig. 8, the present invention includes a plurality of quadrature matched 
filters 60, 62, and 64, which correlate to quadrature matched filters with center frequencies of 
fi, f 2 through f m , respectively, the outputs of which are summed at 66 and compared to a 
threshold at 68. However, unlike conventional post detection diversity combining, or 
averaging, as shown in Fig. 6, the diversity combining 70 occurs prior to detection in the 
present invention. In implementation of the present invention, the received signal r(t) must 
have the transmitter's phase variation removed as fully described hereinbelow. 

Referring to Fig. 9, the validation sequence type diversity combining is nonlinear to 
deal effectively with impulsive noise. Likewise, the differentially coherent combiner must 
contain some nonlinearity to minimize false alarming on impulse noise. Many nonlinear 
filters would work such as median filters, alpha-trimmed filters, and the like. However, to 
maximize processing gain as little data as possible should be discarded. To accomplish this, 
the current implementation of the differentially coherent combiner includes an outlier 
detection algorithm 80 which simply identifies whether all N outputs from the filter are 
reasonably close to one another. If there are a few outliers, they are discarded prior to 
averaging. If there are no outliers, none are discarded. If there are too many outliers (the 
spread of samples is too high), then the whole set of data is discarded as unreliable. 

The outlier detection algorithm 80 can be implemented as follows. First, N samples 
are sorted by magnitude at 81 . If the 3 rd largest sample is much greater than the 4 th largest at 
82, the entire set of samples is discarded as unreliable at 83. Otherwise, if the 2 nd largest 
sample is much greater than the 3 rd largest sample at 84, the two largest samples are 
discarded as unreliable at 85, and the remaining samples are averaged at 86. Otherwise, if the 
1 st largest sample is much greater than the 2 nd sample at 87, the largest sample is discarded as 
unreliable at 88 and the remaining samples are averaged at 86. Otherwise, all of the 
remaining samples are averaged at 86. 
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To implement the inventive "differentially coherent combining" in an EAS receiver, 
the initial conditions on the tag signal due to the transmitter must be constant. A simple way 
to do this is to implement a harmonic transmitter. Instead of having a free running transmit 
local oscillator 6, as shown in Fig. 1, a fixed burst waveform must be transmitted every time. 
5 One way to implement this with a linear transmitter would be to have a transmit waveform 
stored for each frequency: low, nominal, and high. When it is time to send a transmit burst, 
the sequence controller selects which one to send to drive the transmit amplifier. 

When using a switching amplifier, a fixed crystal as the reference to a fractional divider 
to generate the 2-x clock frequency for the switching amplifier can be used. The circuitry 
5 = 10 keeps track of how many cycles are sent out. When the correct number of transmit carrier 
Q cycles are sent out, the transmitter is shut off. Care must be taken in the circuitry so that the 

transmitter starts and ends the same with every transmit burst. 
M When a transmit pulse train of identical bursts is analyzed spectrally, it turns out that 

2 the only signal energy appears at harmonics of the pulse repetition rate, e.g., 90 hertz. Thus, 

W 1 5 even though the transmit energy is centered at 58000 hertz, for example, an infinite pulse 
p train would have zero energy at 58000 hertz. Indeed, the combiner averaging 70, illustrated 

W i n pig. s 5 can b e viewed as a comb filter matched to 90 hertz harmonics. On the other hand, 

m such a combiner will not generally work for a transmitter with a free running oscillator as 

shown in Fig. 1 . In this case, the signal energy does contain 58000 hertz, plus side bands at 
20 integer offsets of 90 hertz from the carrier (due to the amplitude modulation of the 90 hertz 
pulse train). This signal would be heavily attenuated by a 90 hertz comb filter. 

An alternate implementation of differentially coherent combining is to lock the 
receive local oscillator and the transmitter local oscillator in phase and frequency. In this 
way, the carrier phase roll induced by the transmit oscillator would be exactly cancelled by 
25 the phase roll of the receive oscillator. 

The performance of the differentially coherent combining detection scheme of the 
present invention is illustrated as follows. The false alarm probability is again set at 
P fa = 10" 8 . To achieve the same detection probability P de , = 0.968, 14.8 dB SNR is need into 
the noncoherent detector. If N = 4 and receive samples are differentially coherently 
30 combined prior to quadrature detection, we get 3.0*log2 N = 6.0 dB of processing gain. 

Therefore, the raw SNR into the receiver need only be 8.8 dB. This is a 3.2 dB improvement 
over the conventional combining technique. Note the N = 4 is used for convenience of the 
example. In practice N is in the range of 6 to 9. For example, N = 8 gives 9 dB of processing 
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gain. On the other hand, optimum noncoherent combining would give only about 5 dB of 
processing gain. The unanimous vote combiner, which is a suboptimum noncoherent 
combiner, will be even less. In other words, the performance difference becomes greater the 
more diversity is used, the more receive samples are combined. 

It is to be understood that variations and modifications of the present invention can be 
made without departing from the scope of the invention. It is also to be understood that the 
scope of the invention is not to be interpreted as limited to the specific embodiments 
disclosed herein, but only in accordance with the appended claims when read in light of the 
forgoing disclosure. 
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